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Abstract We review some experimental consequences of the presence of super- 
strong magnetic fields of order of the nuclear scale in noncentral heavy-ion colli- 
sions. We present lattice estimates for the strength of the Chiral Magnetic Effect 
(CME) for different quark flavours and argue that the dependence of the anisotropy 
^ • of the distribution of emitted hadrons on their flavor content might be used as an- 

other experimental evidence of the CME. Another possible effect of superstrong 
magnetic field might be the observed abnormal enhancement of dilepton yield. We 
(— i | show that the presence of the magnetic field leads to a specific anisotropy of the 

dilepton emission rate. 



1 Introduction 

It has been realized recently that in heavy-ion collision experiments hadronic mat- 
ter is affected not only by extremely high temperatures and densities, but also by 
superstrong magnetic fields with field strength being comparable to hadron masses 
squared. Such superstrong fields are created due to the relative motion of heavy ions 
themselves, since they carry large charge Z ~ 100 [1]. 

Obviously, the magnetic field is perpendicular to the collision plane, which can 
be reconstructed in experiment from the angular distribution of produced hadrons 
[2]. There is no direct experimental way to measure the absolute value of the field 
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strength, but it can be estimated in some microscopic transport model, such as the 
Ultrarelativistic Quantum Molecular Dynamics model (UrQMD) [ ]. 

Probably the most notable effect which arises due to magnetic fields in heavy-ion 
collisions is the so-called Chiral Magnetic Effect. The essence of the effect is the 
generation of electric current along the direction of the external magnetic field in 
the background of topologically nontrivial gauge field configurations [1, 4]. Such 
generation is not prohibited by ^-invariance, since topological charge density is 
a pseudoscalar field and thus nonzero topological charge explicitly breaks parity. 
However, since QCD is parity-invariant, the net current or the net electric charge 
should vanish when averaged over multiple collision events. Nevertheless, the non- 
trivial effect can still be detected if one considers dispersions of electric current or 
electric charge [5, 6]. Experimentally, the Chiral Magnetic Effect manifests itself 
as the dynamical enhancement of fluctuations of the numbers of charged hadrons 
emitted above and below the reaction plane [7, 2, 8, 9]. 

Another effect, which is closely related to the CME, but has different experi- 
mental signatures, is the anisotropic electric conductivity of hadronic matter in the 
strong magnetic field, which was discovered in lattice simulations [10]. Since the 
conductivity of the hadronic matter is directly related to the lepton emission rate 
[11, 12], such anisotropic conductivity should result in specific anisotropy of the 
dilepton emission rate w.r.t. the reaction plane. This anisotropy should grow with 
the centrality of the collision and with the charges of the colliding ions. This effect 
might also contribute to the observed abnormal dilepton yield in heavy-ion colli- 
sions [ ]. Some theoretical considerations [14] as well as preliminary lattice data 
[15] suggest that at very strong magnetic fields with strength eB > mp\ where m p is 
the mass of the charged p -meson, the anisotropic conductivity might even turn into 
the anisotropic superconductivity in the direction of the field. Unfortunately, such 
extremely large field strengths are hardly reachable with present-day heavy-ion col- 
liders. 

In this paper we give some estimates of the expected experimental signatures of 
superstrong magnetic fields, basing on the lattice data. In Section 2 we consider the 
Chiral Magnetic Effect and argue that its strength should decrease with increasing 
quark mass, which can be used to discriminate between the CME and other pos- 
sible effects which might result in preferential emission of charged hadrons in the 
direction perpendicular to the reaction plane. In Section 3 we consider the dilepton 
emission rate, and estimate the contribution of the induced conductivity to the total 
dilepton yield and dilepton angular distribution in heavy-ion collisions. 



2 Chiral Magnetic Effect 

Chiral Magnetic Effect is usually characterized by the following experimental ob- 
servables, suggested first in Ref. [8]: 
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I N e i N a N b 

a "b = fI E E cos (&« + *Jb) > (!) 

where a,b = ± denotes hadrons with positive or negative charges, respectively, N e 
is the number of events used for data analysis, N a and Nb are the total multiplicities 
of positively/negatively charged particles produced in each event, and 0,- a , 0^ are 
the angles w.r.t. the reaction plane at which the hadrons with labels ; and j are emit- 
ted. Summation in (1) goes over all produced hadrons. In practice, only sufficiently 
energetic particles are considered. 

The main signatures of the CME is the growth of a u b with impact parameter, 

as well as the negativity of a ++ and a and the positivity of a + _ [7, 2, 8, 9]. 

However, it has been pointed out recently that these results can also be explained 
by other effects, such as the influence of nuclear medium on jet formation [16] and 
other in-medium effects [17]. It is therefore important to think about more refined 
experimental tests of the CME. Our main message in this paper is that the depen- 
dence of the charge fluctuations on quark mass can be used to discriminate between 
the CME and other possible phenomena which contribute to the observed asymme- 
try of charge fluctuations. Indeed, the CME emerges due to fluctuations of quark 
chirality [1], which are suppressed when the quark mass is increased. 

The dependence of the observables (1) on the quark mass can be studied if one 
sums separately over charged mesons with different quark content, e.g. ud and du 
(charged pions), us and us (charged kaons) or uc, cu, dc, dc (D-mesons). The ob- 
servables a a j, should then decrease with the meson mass. The dependence of a a b 
on the centrality of the collision should also become weaker. Here we give a rough 
estimate of this effect basing on the results of lattice simulations. 

Let us first note that the observables a a b can be expressed in terms of the dif- 
ferences of multiplicities of charged hadrons emitted above and below the reaction 
plane [1]: 

_ c{A a A b ) 

aab ~ (N a ){N b )> {I) 

where a, b = ±, A± are the differences of the multiplicities of hadrons with positive 
or negative charges above or below the reaction plane, respectively. The factor c 
depends on the hydrodynamical evolution of hadronic matter, and is usually close 
to unity. For multiplicities (N a ) <~ 1000 one can also neglect the initial charge of 
heavy ions Z ~ 100 with a good precision, and assume that (N a ) = (Nb) = N q , 
where N q is the mean multiplicity of the same-charge hadrons per event. 

Lattice results can be compared to experimental data by considering the quantity 

a ++ + a — 2a + _, which can be expressed solely in terms of the difference A Q = 

A + — A- of net charges of hadrons emitted above and below the reaction plane: 

((AQ) 2 ) ((A + -A_) 2 ) 
a ++ + a — -2a+-= 2 = — 2 , (3) 
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In turn, the dispersion of the charge difference ( (A Q) 2 ) can be related to the vacuum 
expectation values of the squared current densities (ft (x)) [5]. The contribution 

of each quark flavor / = u,d,s,c to the total electromagnetic current is (x) = 
qf y^q-t . We do not consider here the third-generation quarks, which are extremely 
rarely produced in heavy-ion collisions. 

The simplest model which allows to express ( (AQ) 2 ) in terms of (j^ 2 (x) ) is the 
model of spherical fireball, which emits positively and negatively charged hadrons 

f 2 

from its surface with intensity proportional to (yj, (x) ). This leads to the following 
relation [5]: 

« ++ + ---^-^^(((;0 2 ) + 2((;{) 2 )), (4) 

where jt (x) and jj_ (x) are, respectively, the currents along the magnetic field and 
perpendicular to it, T is some characteristic collision time, r is the fireball radius 
and p is some typical correlation length for electric charge density in the fireball. In 
our estimates, we take T ~ 0.3 fm (this is a typical decay time for the magnetic field 
in heavy-ion collisions [3]), p ~ 0.2 fm, which are reasonable parameters for, say, 
gold-gold collisions at 60GeV/nuclon. We also assume that the fireball is a sphere 
with radius r = R — b/2 within the overlapping region between the two heavy ions of 
radius R which collide at impact parameter b. The net multiplicity N q and the impact 
parameter b as the functions of collision centrality can be found in Table 1 in [1] . For 
simplicity, we assume that ( (j^)^ 2 } are approximately constant on the surface of 
the fireball. Note also, that in order to exclude the effects related to the dependence 
of multiplicities of strange and charmed mesons on the collision centrality, which 
might be different from that of light mesons, we normalize the charge of emitted 
hadrons by the square of the total multiplicity of all hadrons. 

Several technical remarks are in order. First, we assume that the matter within the 
fireball is in the state of thermal equilibrium, and thus the expectation values (j^ 2 ) 
can be calculated from gauge theory in Euclidean space. We also assume that the 



Fig. 1 Schematic view of 
the collision geometry. The 
fireball is the hatched region 
of volume V ~ ^ (R - b/2) 3 
within the intersection of two 
heavy ions of radius R each. 
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magnetic field is uniform and nearly time-independent. Of course, these are rather 
rough approximations, but we are aiming here at qualitative rather than quantitative 
estimates. We have calculated the currents in SI/ (2) lattice gauge theory with back- 
ground magnetic field both in the confinement phase, neglecting the contribution 
from the virtual quark loops (quenched approximation). A comparison with SU (3) 
gauge theory suggests that this is a reasonable approximation [ 1 8] . A more detailed 
study has shown also that in the deconfinement phase the dispersions of local cur- 
rent densities are practically independent of the magnetic field [5], thus we do not 
consider here this case. The expectation values (jjj 2 } contain also the ultraviolet di- 
vergent part, which we have removed by subtracting the corresponding expectation 
values at zero temperature and zero magnetic field. 

The masses of the valence quarks took the values m q = 50MeV (for u and d 
quarks), m q = WQMeV (for s-quark) and m q = IGeV (for c-quark). The lowest 
value of the quark mass is dictated by the numerical stability of our algorithm. Thus 
our calculations of ( j^ 2 } at m q = 50MeV should be considered as only the lower 

bound for (j^ 2 ) at realistic masses of u and d quarks. 

On Fig. 2 we compare the quantity a ++ + a — 2a + _ calculated for the experi- 
mental data by the STAR collaboration [ ] with the estimate (4) based on the results 
of lattice simulations at different quark masses and different temperatures. In order 
to match the value of the magnetic field strength in experiment and in simulations, 
we take the rough estimate from eq. (A. 12) in [1]: 

eB~ (0.1£/fl)GeV 2 , (5) 



10 
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Fig. 2 Comparison of the quantity a ++ + a — 2a^ for the experimental data by the STAR 

collaboration [2] with the estimates (4) based on the results of lattice simulations at different quark 
masses. 
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where b is the impact parameter. This rough fit also agrees by the order of magnitude 
with the results of more sophisticated calculations of the magnetic field within the 
UrQMD model [3]. 

One can see that the best agreement with the STAR data is obtained at the small- 
est quark mass. The combination a ++ +a — 2a H quickly decreases as the quark 

mass is increased - approximately by a factor of 5 as the quark mass changes from 
SQMeV to 1 IQMeV, and by a factor of almost 20 as the quark mass further increases 
to 1 GeV . Since all observables a a b are typically of the same order, one can expect 
that each such observable will also decrease with the quark mass. This dependence 
of asymmetry of angular distributions of mesons of different flavors on their mass 
can be used to discriminate between the CME and other phenomena which might 
cause such asymmetry. 

The result may be compared with the perturbative analog of CME [ 1 9] resulting 
from straightforward generalization of Heisenberg-Euler Lagrangian depending on 
quark mass as m~ 4 



The correspondence of this perturbative and Abelian effect to CME is manifested 
by the substitution 



As a result, the perturbative mechanism may become essential when quark mass is 
exceeding the inverse correlation length of topological charge density. One may ex- 
pect that the transition point from non-perturbative to perturbative mass dependence 
is not too far from the strange quark mass, like it happen for vacuum quark conden- 
sates and strangeness polarization in nucleons (see [19] and Ref. therein) so that the 
non perturbative lattice results are applicable for experimentally important case of 
strangeness separation. 



3 Induced conductivity and abnormal dilepton yield 

Another phenomenon which might be caused by superstrong magnetic fields act- 
ing on the hadronic matter is the induced anisotropic conductivity along the mag- 
netic field [10]. While the Chiral Magnetic Effect is related to the local fluctuations 
of current density, the induced conductivity reflects the fact that these fluctuations 
also have long-range correlation in time. Indeed, by virtue of the Green-Kubo rela- 
tions the conductivity is related to the zero-frequency limit of the spectral function 
Pfiv (w) which corresponds to the correlator ( & (x) j v (y) ) in Minkowski space 




(GG) 



(6) 




(7) 



[20]: 
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dw coshw (t— 



Jd 3 x{^j^ (0,0), v (x,T)> = j 



In sinh(^) 



PflV (w) 



(8) 



o 



This spectral function can also be extracted from the results of lattice simulations 
using the so-called Maximal Entropy Method [21, 22]. 

The spectral function p^ v (w) determines also the dilepton emission rate from 
either cold or hot hadronic matter [23, 11]: 



is the dilepton tensor (T7^ v is the Minkowski metric), m is the lepton mass. Thus the 
low-momentum limit of p^ v (w) is related, on the one hand, to the emission rate of 
soft dileptons, and, on the other hand, to the conductivity of hadronic matter. 

The enhancement of the conductivity due to the magnetic field should thus lead 
to the enhancement of the dilepton emission rate. This might provide a viable ex- 
planation of the abnormal soft dilepton yield observed in heavy-ion collisions [13]. 
Moreover, the anisotropy of the conductivity should lead to specific correlations be- 
tween the momenta of the dileptons and the direction of the magnetic field (in other 
words, with the orientation of the reaction plane). 

In a similar way, one can also consider the correlators of charged vector currents 

j{? = {q^YiiQ^ )> where f =/= f are the flavour indices. The correlators of charged 
vector currents have been studied recently in [15] in order to test the conjectured 
p -meson condensation in superstrong magnetic fields [14]. It was found that as the 
magnetic field is switched on, this correlator decays slower, which corresponds to 
the the increase of the corresponding spectral function at lower frequencies. Since 
charged currents are coupled to leptons and neutrinos, this effect might be man- 
ifested in the enhanced lepton and neutrino yield. Excitations over the p -meson 
condensate might show up then in the usual decay channel p ± — s- n ^. 

Lattice simulations show that the electric conductivity is nonzero only in the 
direction of the magnetic field and depends linearly on qB [10, 24], For suffi- 
ciently small momenta p\ and pi (and thus for small q) by virtue of the Green- 
Kubo theorem one has Pij(q) ~~ ^ijl/T ~ BjBjq/ (\B\T) [20, 21]. Let us also 
neglect the lepton masses and go to the rest frame of the dilepton pair, where 
Pi = — P2 = pn, q = (2p,0) and the spatial components of the dilepton tensor are 
L' 7 = p 2 (8ij — n,n,). The dilepton emission rate is therefore proportional to 




(9) 



where p\ and p2 are the momenta of the leptons, q = p\+ pi, 





(10) 
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where 9 is the angle between the spatial momentum of the outgoing leptons and the 
magnetic field. Therefore, there should be more soft dileptons emitted perpendicular 
to the magnetic field than parallel to it. As a result, they are to large extend hidden 
inside the hadrons in the scattering plane which should lead to the difficulty in their 
experimental observation. Similar predictions can be also made for the angular dis- 
tributions of 7r ± 7T° pairs in the case of charged meson decays (see above). 

In order to estimate the effect of the magnetic field on the total dilepton yield 
in heavy-ion collisions, we normalize the conductivity induced by the magnetic 
field to the conductivity at zero magnetic field and at the temperature close to the 
deconfinement phase transition, T = \.\2T C . By virtue of (9), the ratio of these 
conductivities should be equal to the ratio of dilepton emission rates in the low- 
momentum region. The relevant lattice data is summarized in [10, 24]. In these 
works it was found that in the deconfinement phase the conductivity is practically 
independent of the magnetic field. Therefore here we will try to estimate possible 
contribution of the magnetic field to the dilepton emission rate from hadronic matter 
in the confinement phase. From the data presented in [10, 24], we estimate the ratio 
of the induced conductivity a (B,T < T c )to the conductivity of quark-gluon plasma 
C7(B = 0,r = 1.12r c ) as: 



We now use the estimate (5) for the magnetic field strength and take into account 
that the ratio of conductivities should be equal to the ratio of dilepton emission rates. 
We thus obtain for the contribution of the magnetic field to the dilepton emission 
rate: 



where the factor 1/3 appears after averaging the expression (10) over the angle 9 . 

We conclude that at large impact parameters (b ~ 2R) the dilepton yield can in- 
crease by up to 20 — 30% due to the influence of the magnetic field. This factor 
should be essentially reduced because the magnetic-induced dileptons are hidden 
inside the scattering plane. The observed abnormal dilepton yield is, however, max- 
imal for central collisions (where it reaches several hundred percent as compared 
to the hadron resonance model) and decreases as the impact parameter grows [13]. 
Such behavior might be caused by several factors, such as the change of the temper- 
ature within the fireball or the change of the fireball volume with impact parameter. 
A proper investigation of such factors is out of the scope of this paper, and cannot 
be undertaken using the methods of lattice gauge theory. 



<t(B,T < T c ) /a (B = 0, T = 1 .12r e ) 



eB 



(0.5 GeV) 2 




(11) 



Lattice studies of magnetic phenomena in heavy-ion collisions 

4 Conclusions 



9 



In this paper we have summarized the main experimental signatures of the effects 
caused by superstrong magnetic fields in heavy ion collisions, namely, the Chiral 
Magnetic Effect and the abnormal dilepton yield. 

The Chiral Magnetic Effect [ I ] results in preferential emission of charged hadrons 
in the direction perpendicular to the reaction plane. The origin of the Chiral Mag- 
netic Effect is the fluctuations of chirality, which are suppressed as the quark mass 
grows. Thus this asymmetry in angular distributions of charged hadrons, character- 
ized by the coefficients a a \ (1), should be strongly suppressed for strange or charmed 
hadrons. 

The abnormal dilepton yield with specific angular dependence (10) is the con- 
sequence of electric conductivity of the hadronic matter induced by the magnetic 
field. In this case more dileptons are emitted in the direction perpendicular to the 
magnetic field. Let us also note that since according to our lattice data the magnetic 
field influences the conductivity only in the confinement phase [10], the significant 
change of dilepton yield in noncentral heavy-ion collisions might be a signature of 
the confinement-deconfinement phase transition. 

More generally, lattice data suggests that the influence of the magnetic field on 
the properties of hadronic matter is stronger in the confinement phase. Thus heavy- 
ion collision experiments on colliders with lower beam energy but with larger lumi- 
nosity (such as FAIR in Darmstadt, Germany or NICA in Dubna, Russia) might be 
more advantageous for studying magnetic phenomena. 
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